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Abstract 

A possibility that in the FPU problem the critical energy for chaos goes to 

zero with the increase of the number of particles in the chain is discussed. 

The distribution for long linear waves in this regime is found and an estimate 

for new border of transition to energy equipartition is given. 
PACS numbers: 05.45, 63.70 



Typeset using REVTpjX 



Starting from 1955 the Fermi-Pasta-Ulam (FPU) problem JlJ initiated numerous re- 
searches and became one of the corn stone in the modern statistical mechanics . The 
absence of energy equipartition in the system of coupled nonlinear oscillators observed nu- 
merically in ]T| pushed forward the investigations of chaos as well as the analysis of com- 
pletely integrable systems (see || and refs. therein). 

The first explanation of the striking result 0] was proposed by Chirikov and Izrailev 
Q on the basis of the Chirikov criteria of overlapping resonances [||]. According to [§]] it 
is necessary to exceed some critical energy value to obtain overlapping of the resonances, 
chaos and energy equipartition over linear modes. According to |Q in the case of low mode 
excitation (nonlinear sound waves) the critical energy increases with the number of oscillators 
in the chain (or energy per oscillator is constant). Below this energy it was argued that the 
resonances are not overlapped and the motion is close to integrable one. Since some of the 
initial conditions in |l| were below this border the energy equipartition was absent ||]. The 
results of ||] were confirmed in the series of analytical and numerical researches || where the 
authors also analyzed the dependence of the Lyapunov exponents on the energy. However, 
these researches showed that the relaxation to an equilibrium distribution could be very long 
at small energies that makes difficult to study the transition from global chaos to intergable 
case. 

In this paper the condition of resonance overlapping for long (sound) waves in the limit 
of small energy is analyzed. For long waves the dispersion law is very close to linear. Due 
to that for the system with finite but large number of oscillators N there are some terms in 
the nonlinear part of the Hamiltonian which are in the resonance even for very low energies. 
Such resonances being not considered in |4j] give a sharp decrease of the chaos border in 
energy which goes to zero with the increase of the number of particles in the lattice. In this 
sense the long wave chaos can exist for arbitrary small nonlinearity. The physical reason 
of this unusual phenomenon is connected with the linearity of unperturbed system. Due to 
that for the sound dispersion law which is typical for long waves the KAM theorem cannot 
be applied and chaos can appear for arbitrarily small nonlinear interaction. Such kind of 
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phenomenon have been already studied in different dynamical systems with few degrees of 
freedom |7]-|9| • 111 such a case the dynamics can be described by a renormalized Hamiltonian 
independent on the strength of nonlinear interaction. In particular the measure of chaotic 
component remains unchanged with the decrease of nonlinearity. In the FPU-problem a 
deviation of the dispersion law from the linear one gives rise to a critical chaos border which 
is, however, extremely low and decreases with the number of particles in the chain. 

We start our analysis from the a-FPU problem with cubic nonlinearity in the Hamilto- 
nian: 

1 N a N 

iJ= aEb» 2 + ~ X nf] + ( X n+1 ~ X nf (1) 

where the first sum gives the Hamiltonian H of the linear waves and the second sum 
represents the interaction Hi nt . The boundary conditions are fixed as xq = 0;xn+i = 0. 
The eigenmodes (Q k , P k ) of H Q are connected with the coordinates x n ,p k by the equations 
x n = v /2/(ivTl)Efe<5jfcSin(g fc n), Pn = p/(N + 1) £ fe P k sin(q k n) with q k = irk/(N + l), 
1 < k < N 0. In this representation H can be written as H = J2k (Pk 2 + ^fc 2 Qfc 2 )/2 = 
J2k u klk with the eigenfrequencies u k = 2sin(g fc /2). The action-angle variables (Ik, Ok) are 
connected with (Pk,Qk) in the standard way 0. 

It is convenient to write the total Hamiltonian in the action-angle variables (Ik, 6k) of 
the linear problem. Taking into account that the nonlinear coupling is small we can keep 
in H int only the resonant terms corresponding to the resonant 3-waves interaction. For 
long waves this condition corresponds to ^3 = ^2 + ^1- All other terms can be eliminated 
by averaging over fast oscillations with frequencies iOk- After this procedure we obtain the 
averaged Hamiltonian: 

H = ^uJkh + y (^k^^kjkj^hs) 112 cos(8 k3 - 6k 2 - 9 hl )S k3tkl+k2 (2) 

k 2 viV + 1 kiMM 

which can be written as H = Hq + Hi nt . Here the bar marks the averaging over fast 
oscillations with uj k > oo\ = ir/(N + 1). The term fast means that uo\ ^> 5u where 5u is 
the typical nonlinear frequency 5cu ~ dH int /d9 ~ a(E /N) 1 ^ 2 tUk- For few low modes excited 
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around a given fc-value we obtain 5uj ~ a (Eo/Nf^k/N where E is initial energy. Following 
the way of |4[], where the /3-FPU model with quartic nonlinearity had been studied, we can 
find the chaos border from the condition of the overlapping resonances 5uj ~ Auj where 
Au w uj\ ~ Ti /N is the distance between the main resonances in (1). According to this 
condition the global chaos appears for at = aE Q 1 ^ 2 > olchi ~ v^V /k. 

The Hamiltonian H has additional integral of motion Eg = TvY^kklk/ '(N + 1) ~ E Q . 
For long sound waves (k << N) we can use approximate expression for the dispersion law 
^fc = Qk — 9fc/24 in H while in the term with H int it is sufficient to use oj k = %■ By using the 
new resonant phases <p k = Ok — Qkt we can transform (2) to the new resonant Hamiltonian: 

M M M-fei 

H R = _cr H k z l k + 2/i {kik2k k2+kl I k Jk 2 h 2 +k 1 ) l/ ' 2 cos(0 fe2+fcl - 4> k2 - (fikj (3) 

k=l fei=l k 2 =l 

where a = 7r 3 /(24(iV + l) 3 ), \i = 7r 3//2 a/(4(iV + l) 2 ) and M is the maximal number of 



harmonics. It is convenient to introduce the new dimensionless time r = fityEs(N + 
in which the dynamics is described by the renormalized resonant Hamiltonian 

Hrn = —vJ2 k =i k^Jk , . 

(4) 

+2E"=i EK=i fcl (^ife(fc 2 + fci)J fel J fc2 J fc2+fcl ) 1 / 2 cos(0 fe2+fcl -<p k2 -<f> kl ) 



with one dimensionless parameter 



(5) 



Hy/(N + 1)E S Ga^(N +1) 3 / 2 ' 
The new actions J k are connected with the old ones by the relation I k = EgJk(N+ 1)/tt. 
They are now normalized by the condition Y^f=i kJ k = 1. 

Let us analyze now the dynamics of the system (4). If initially only few modes are excited 
around a fc-value then the distance between the resonances is Aoo w vk 2, while the width of 
the resonance is 8u ~ k z ^ 2 J k 1 ^ 2 ~ k. From these estimates it is clear that the resonances 
are overlapped for v < v cr ~ 1/fc 2 and then chaos arises. In the original variables this 
means that the chaos border is given by 

/ — k 2 

a\JE s > a s ^ or v<l/k 2 (6) 
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This border, which takes into account the degenerate sound resonances k 3 = k 2 + hi, 
decreases with the growth of N and is N 2 times below the border of global chaos a C Hi- In 
the case of excitation of low modes with k ~ 1 the critical energy above which the motion 
is chaotic is E c ~ l/(a 2 N 3 ). Therefore, chaos arises at zero temperature T = E /N. 

For a better understanding of the properties of the system (4) a numerical investigations 
of its dynamical motion was carried out. The initial conditions were usually fixed as three 
excited modes with J\ = J 2 = J 3 = 1/6 and different phases 0. The calculations of the 
maximal Lyapunov exponent shows that above the border (6) the motion is characterized by 
the positive exponent \rn that indeed demonstrates the existence of chaos in this regime. 
Above the border the maximal Lyapunov exponent is zero (except exponentially narrow 
chaotic layers). A typical example of the dependence of \rn on renormalized time r is 
presented in Fig. . The energy distribution E k = k.J k over linear modes is shown in Fig. . 
To suppress the fluctuations the values of E k were averaged over time r in the time interval 
[1000-2000]. Below the chaos border (6) the number of excited modes remains the same as 
for the initial distribution. On the contrary above this border the energy is distributed over 
some finite width Ak which is much larger than the initial width. For the high values of 
k ^> Ak the distribution decays in an exponential way. In the whole interval of k the energy 
distribution Ek can be fitted by the effective distribution: 

A 

fk = Z(exp(fc/Z - 7 ) + 1) (7) 
where the length / determines the effective number of excited modes, 7 is some constant 
which mainly effects the shape of the distribution for small k and A is determined by 7 via 
the normalization condition J2Ek ~ / fkdk = 1. For the case of Fig. the optimal value is 
7 = 2.65. It is interesting to note that the fitting (7) quite well describes the distribution 
E k in the large interval 0.03 < v < 0.0005 with the same 7 and different /. This fact is 
demonstrated in Fig. where 6 distributions are superimposed in the rescaled variables lE k 
and k/l. The fitting (7) allows to determine the dependence of the length / on v. This 
dependence is presented in Fig. and is approximately given by the equation I = 0.42/y^. 
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The same functional dependence on v takes place for the quantity 1/ E\ which characterizes 
the width of the distribution Ak for small k. The existence of the same scaling on v for I 
and 1/ Ei confirms once more that the distribution E\~ has only one scaling parameter I. 

The obtained scaling of I from v can be understood on the following grounds. The 
nonlinear resonance width in (4) is 5u ~ dH RN /dJ k ~ k^y^k 1 / 2 with k ~ I. The last 
term k 1 ^ 2 gives the result of summation over k terms with random phases contributing in 
5uj. A typical distance between the resonances is Au; ~ vk? . The number of excited modes 
is determined by the chaos border given by the resonance overlapping: 5u > Au. According 
to this estimate the number of excited modes is Ak ~ I ~ l/v^ that is in agreement with 
the numerical dependence from Fig. and the previous estimate (6). Using the expression 
for v we can find the effective number of excited linear modes expressed via the original 
variables: 

AA; ~ / ~ (a 2 E N 3 f' A (8) 

From this expression it follows that for fixed a and E the number of excited modes is 
quite large but still Ak/N < I. 

In the same way we can obtain estimate for the maximal Lyapunov exponent \rn in 
the renormalized Hamiltonian (4). Indeed, \rn ~ Su ~ k 3 ^ 2 y/k.Jk with k ~ I and \rn ~ 
k ~ Using the relation between the time t for the original system (1) and the time 

r in the renormalized Hamiltonian (4) we obtain the estimate for the maximal Lyapunov 
exponent A in the system (1): 

A(N + l) 3 / 2 ~ [ ' 

The numerical data for the dependence of \rn on v are presented on the Fig. . Un- 
fortunately, in the given interval of v the variation of A^at is not quite monotonic and 
further numerical investigations are required for verification of the theoretical dependence 
Ar/v ~ ( see t ne discussion below). Let us mention that the sign of v in (4) is not 

important and the results are qualitatively the same for v < when the absolute value of v 
should be used in the estimates. 
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The comparison of A with the distance between main resonances Au; shows that for 
sufficiently large N the nonlinear resonance width 8u ~ A becomes larger than Au; ~ 1/N. 
The condition A > Au shows that the main resonances in (1) will be overlapped for 



Above this border the nonresonant terms neglacted in the derivation of (2) give the 
overlapping of the main resonances and for aEo > 1/N 1 ^ 3 approximate equipartition over 
all linear modes modes can be expected. So, in the limit of large N the equipartition can 
appear at zero energy and zero temperature. The time required to reach the equipartition 
is inversely proportional to A. 

It is interesting to note that some conditions of considered usually as integrable 
(Fig.l in |J) have v m 0.13. Direct computation in (4) for this v value with corresponding 
initial conditions gives, however, Xrn = 0. This puts the question about a more exact 
determination of the border of chaos \v CT \. 

Let us now briefly discuss the properties of chaos in the /3-FPU model with quartic 
interaction H int = PY^ n ( x n+i — ^n) 4 /4. As in the a-case we should keep only the resonant 
terms for 4 waves with k\ + A; 2 = k 3 + /c 4 . The resonance nonlinear width can be estimated in 
the same way as in j|j2| 5u ~ {3Eo(j}./N. The overlapping of the main resonances happens for 
(3E > N/k However, for the resonant Hamiltonian only the deviation of Uk from the 

sound law nk/N is important so that the distance between the resonances can be estimated 
as Au ~ k 3 /N 3 [|nj. This gives the border of slow chaos (3E > k 2 /N which is much below 
the standard border @,@]. Above this border the number of excited low linear modes is 
k m Ak ~ ^/(3EqN and the maximal Lyapunov exponent is A ~ (PEq/N) 3 ^ 2 ~ 8u. The 
overlapping of the main resonances takes place for A > 1/N or (3Eq > N 1 ^ 3 . Above this 
border all linear modes are excited leading to energy equipartition. In a difference from the 
a-model this border grows with N but the critical temperature T = Eq/N still goes to zero. 

The above theoretical estimates were based on the comparison of the splitting between 
linear modes and nonlinear spread width. As in the case of the Chirikov criteria such 




(10) 
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approach cannot exclude a possibility that the system under investigation is completely 
integrable or is very close to some of them. This point is very crucial for the a-FPU problem 
since at low energy it is very close to the Toda lattice (see 0). Due to that generally we 
should expect that contrary to the above estimates and numericaly data the dynamics of 
a-FPU problem will be integrable. To understand this apparent contradiction with the 
numerical data additional simulations had been carried out. Namely, the total number of 
harmonics M had been increased up to M = 120 for the parameters of Fig. la. While the 
simulations become very heavy in such a case they give approximately In r/r decay of Xrn 
up to A ~= 0.02 at maximally reached r = 400. This indicates that in a real system with 
very large M the Lyapunov exponent will be zero. At the same time such change of M did 
not affected the averaged energy distribution (see Fig. 2). For a better check of this point a 
number of numerical simulations with the original Hamiltonian (1) had been done with TV 
up to 151 and the initial conditions corresponding to the Fig. la with fixed v = 0.01. For 
N = 61 the renormalized Lyapunov exponent (see (9)) was stabilized around Xrn ~ 0.13 
(the time t max in the simulations was t max ~ 9 x 10 5 ); for N = 101 the exponent was also 
stabilized around Xrn ~ 0.065(t maa; ~ 9 x 10 6 ). However, in both these cases the averaged 
energy distribution had significant increase at high modes in a difference from Fig.2. For 
N = 151 during all t < t max ~ 6 x 10 6 the value of \rn was decreasing as In r/r reaching 
Xrn ~ 0.04 at t max . At the same time the averaged distribution was practically the same 
as on the Fig.2 (see Fig. 5). These additional data show that in the low energy limit the 
dynamics of a-model is not chaotic (Xrn = 0) as it can be expected from the comparison 
with the Toda lattice. Inspite of that the energy distribution (see (7)) is correctly given by 
the above estimates derived from the renormalized Hamiltonian. The reason due to which 
the renormalized dynamics is so sensitive to the maximal value of M is still not quite clear. 
It is possible that the important effects of coupling to high modes can be understood from 



nonlinear wave equation in the continuous limit (see |Tl[). Very recently, the properties of 



the Lyapunov exponent in the system (1) with N up to 128 were studied in |T2| . 

The situation for the /3-model can be more interesting. Indeed, apparently this model is 
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not close to any integrable system and the above renormalization approach and estimates 
should give correct chaos border. The picture of low energy chaos developped here is quali- 
tatively close to that one in ||13| . However, additional investigations of this regime are still 
highly desirable. They should clarify some uncertainties in the estimate of Au> (see [pHH). 
Also the question of coupling to high modes can play a very crucial role [ini . 

I would like to thank J.Bellissard and P.L'Eplattenier for stimulating discussions. The 
above results had obtained at the begining of 1994. I greatly appreciate the discussions and 
critisism of B.V.Chirikov being in Toulouse at that time. 
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FIGURES 



Fig. 1. Maximal Lyapunov exponent \ R n m (4) as a function of time r: a)full line: 
v = 0.01, H RN = 0.544, \ RN > 0; b)dashed line: v = 1., H RN = —5.45, \rn — > (values 
of \ R n are multiplied by 5). For all Figs. 1-4 only three modes were initially excited with 
j 1 = j 2 = j 3 = i/6. 

Fig. 2. Averaged energy distribution E k = kJ k over linear modes k for the cases of 
Fig.l: a) full circles; b) open circles. Full line gives the fitting distribution (7) with / = 4.05; 
7 = 2.65, A = 0.3678. 

Fig. 3. Normalized energy distribution lE k as function of k/l for the six different values 
of v from Fig. 4 (points). The full line gives the fitting (7) with 7 and A from Fig. 2. 

Fig. 4. Dependence on v for: length / obtained from distributions of Figs. 3 (points); 
average energy of first mode E 1 (squares); \ RN (open circles). The straight line shows the 
theoretical dependence / ~ 

Fig. 5. Same as Fig.2 obtained from the Hamiltonian (1) (see text). 
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